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A new version of a Hartree-Fock program is presented that includes ex- 
tensions for partially filled /-subshells. The program allows the calculation 
J^i | of term dependent Hartree-Fock orbitals and energies in LS coupling for 

^^ ' configurations with no more than two open subshells, including /-subshells. 
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PROGRAM SUMMARY 

Title of program: HF96 

Catalogue number: ADDZ 

Program obtainable from: CPC 
Program Library, Queen's 
University of Belfast, N. Ireland (see 
application form in this issue) 

Licensing provisions: none 

Computer: SUN Sparcstation 1 

Installation: Vanderbilt University, 
Nashville, TN 37235, USA 

Operating System under which the 
program is executed: Sun UNIX 

Programming language used: 
FORTRAN 77 

Memory required to execute with 
typical data: 449K Bytes 

No. of bits in a word: 32 

Peripherals used: terminal, disk 

No. of lines in combined program 
and test deck: 4703 

CPC program Library subprograms 
used: none 

Keywords: atomic structure, 
complex atom, Hartree-Fock 
method, wave functions, bound 
state, energy level, LS coupling. 



Nature of physical problem 
This program determines 
non-relativistic, numerical radial 
functions and computes energy 
levels in the Hartree-Fock 
approximation . 

Method of solution 
The new version of the program 
retains the design and structure of 
the previous one [1]. The same 
self-consistent field method is used 
to solve the differential equation [2]. 
Tables from Ref.[3] that define 
deviations of term energies from the 
average energy of the configuration 
E av have been extended to include 
the partially filled /-subshells. 

Restriction on the complexity of the 
program 

Only one or two open subshells are 
allowed. Any number of s, p, d and 
/ electrons are allowed in the case 
of one open subshell. In the two 
open subshell case, the program can 
calculate the term dependent energy 
expression only for the following 
configurations: 

1. (ns)(n'l) N ,where 1 = 0,1, 2, 3, 

2. (np) N (n'l) , where 
Z = 0,1,2,3,... , 

3. (nf)(n'd). 

Typical running time 

The time required to execute the 

test case on a SUN Sparcstation 1 



was 64.5 seconds. 
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LONG WRITE-UP 

1 Introduction 

Recently there has been a tremendous growth of interest in lanthanide and 
actinide chemistry and physics, yet relatively few calculations have been 
performed for these sytems. Though accurate results will require a fully 
relativistic treatment in the jj-coupling scheme, that also accounts for the 
nuclear effects on the energy level structure as is done in GRASP92 [1], for 
example, much can be learned about quantitative behaviour of atomic prop- 
erties from simpler non-relativistic calculations. For example, Tatewaki et 
al [2] recently found that the non-relativistic Hartree-Fock ionization poten- 
tials of singly ionized lanthanides follow the observed trends, 

The numerical Hartree-Fock program, HF86, published earlier [3] could 
only treat configurations with open /-subshells containing more than one 
electron, in the average-energy-of-CONFIGURATION approximation [4]. 
In the present version, HF96, this restriction is removed, though like the ear- 
lier version, it is only able to compute the term-dependent orbitals when an 
expression for the energy can be determined from relatively simple formulas. 
More complex configuration states will still need the full power of the angu- 
lar momentum programs part of the MCHF Atomic Structure Package [5] , 
extended also for open /-subshells. 

2 Program organization 

The new version of the HF program has the same design as the previous 
one [3]. It gets an expression for the average energy from the occupation 
numbers for the subshells and then, corrects the average energy expression 
with a deviation so that 

E{LS) = E av + AE(LS). 

The deviation, AE(LS), has the property that the expression is the same 
for particles and holes. Thus one needs to be concerned only with up to 
half-filled subshells and their deviations. There are simple rules for other 
cases (for details see [4]). 

The HF program [3] is extended to include the cases: 



1- (nf) 



N 



2. (nf) N {n's), 

3. (nf)(n'd), 

Input data for the new version of HF is the same as before. One needs 
only to take into account that the classification of terms of the / subshell 
is more complicated than for s, p, d subshells. It is related to the fact that 
for the classification of /-subshell states the characteristics S' (multiplicity), 
L (total orbital momentum), and v (seniority) are not enough. We use a 
classification of terms similar to that used by Nielson and Koster [7] . Namely 
for a particular subshell instead of v, when the term is not unique, a single 
character "number", Nr, is used. Its value is found in Table 1 where all 
terms for /-subshells are presented. In most cases, Nr appears to be a 
digit, but since it is a single character, the single letter A is used instead of 
the number 10. 

In many instances, the term can be uniquely determined and neither 
seniority v nor Nr are need. However, when more than one coupling is 
possible, additional information is requested. Consider a calculation for the 
A F term of the configuration 4f (4)6s(l) for which the coupling of 4f(4) 
could be either 3 F or 5 F. Table 1 shows there is one 5 F term and four 3 F 
terms. In cases of such ambiguity, the user will need to provide additional 
information for the term and Nr of the "parent", or more precisely the l**n 
subshell. This is illustrated in the Test Run Output. 

3 Example 

Our example shows the input to the interactive Hartree Fock program and 
some of the intermediate output. The calculation is for the total energy 
of the level ls 2 2s 2 2p 6 3s 2 3p 6 4s 2 4p 6 4d 10 5s 2 5p 6 4/ 4 3 ^ =(211)c/=(21) 6s 1 4 F of 

Nd + . Initially, calculations were requested for 4 F. Then the program deter- 
mined an Ambiguous l**n parent term and requested information about 
the term and Nr for 4f(4). A calculation was desired for the 3 F parent with 
Nr =3. This information was entered as a three-character string, 3F3. With 
this information, the energy expression could be determined. 



Table 1: Classification of /- subshell states. 
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(20) 


5 


2 D 


3 


(210) 


(20) 


1 


subshell f 7 






4 G 


7 


(220) 


(21) 


6 


2 D 


3 


(210) 


(21) 


2 


S S 


7 


(000) 


(00) 





4 G 


7 


(220) 


(22) 


7 


2 D 


5 


(221) 


(20) 


3 


6p 


5 


(110) 


(11) 





4 tf 


5 


(211) 


(11) 


1 


2 D 


5 


(221) 


(21) 


4 


e D 


7 


(200) 


(20) 





4 H 


5 


(221) 


(21) 


2 


2 D 


5 


(221) 


(31) 


5 


6 F 


5 


(110) 


(10) 





A H 


5 


(221) 


(30) 


3 


2 D 


7 


(222) 


(20) 


6 


6 G 


7 


(200) 


(20) 





A H 


7 


(220) 


(21) 


4 


2 D 


7 


(222) 


(40) 


7 


6 H 


5 


(110) 


(11) 





A H 


7 


(220) 


(22) 


5 


2 F 


1 


(100) 


(10) 


1 


H 


7 


(200) 


(20) 





4 / 


3 


(111) 


(20) 


1 


2 F 


3 


(210) 


(21) 


2 


l S 


3 


(111) 


(00) 


1 


4 / 


5 


(211) 


(20) 


2 


2 F 


5 


(221) 


(10) 


3 


A S 


7 


(220) 


(22) 


2 


4 / 


5 


(211) 


(30) 


3 


2 F 


5 


(221) 


(21) 


4 



LS 


V 


W 


U 


Nr 


LS 


V 


W 


U 


Nr 


2 F 


5 


(221) 


(30) 


5 


2 I 


5 


(221) 


(30) 


3 


2 F 


5 


(221) 


(31) 


6 


2 I 


5 


(221) 


(31) 


4 


2 F 


5 


(221) 


(31) 


7 


2 I 


5 


(221) 


(31) 


5 


2 F 


7 


(222) 


(10) 


8 


2 I 


7 


(222) 


(20) 


6 


2 F 


7 


(222) 


(30) 


9 


2 I 


7 


(222) 


(30) 


7 


2 F 


7 


(222) 


(40) 


A 


2 I 


7 


(222) 


(40) 


8 


2 G 


3 


(210) 


(20) 


1 


2 I 


7 


(222) 


(40) 


9 


2 G 


6 


(210) 


(21) 


2 


2 K 


3 


(210) 


(21) 


1 


2 G 


5 


(221) 


(20) 


3 


2 K 


5 


(221) 


(21) 


2 


2 G 


5 


(221) 


(21) 


4 


2 K 


5 


(221) 


(30) 


3 


2 G 


5 


(221) 


(30) 


5 


2 K 


5 


(221) 


(31) 


4 


2 G 


5 


(221) 


(31) 


6 


2 K 


5 


(221) 


(31) 


5 


2 G 


7 


(222) 


(20) 


7 


2 K 


7 


(222) 


(30) 


6 


2 G 


7 


(222) 


(30) 


8 


2 K 


7 


(222) 


(40) 


7 


2 G 


7 


(222) 


(40) 


9 


2 L 


3 


(210) 


(21) 


1 


2 G 


7 


(222) 


(40) 


A 


2 L 


5 


(221) 


(21) 


2 


2 H 


3 


(210) 


(11) 


1 


2 L 


5 


(221) 


(31) 


3 


2 H 


3 


(210) 


(21) 


2 


2 L 


7 


(222) 


(40) 


4 


2 H 


5 


(221) 


(11) 


3 


2 L 


7 


(222) 


(40) 


5 


2 H 


5 


(221) 


(21) 


4 


2 M 


5 


(221) 


(30) 


1 


2 H 


5 


(211) 


(30) 


5 


2 M 


5 


(221) 


(31) 


2 


2 H 


5 


(221) 


(31) 


6 


2 M 


7 


(222) 


(30) 


3 


2 H 


5 


(221) 


(31) 


7 


2 M 


7 


(222) 


(40) 


4 


2 H 


7 


(222) 


(30) 


8 


2 N 


5 


(221) 


(31) 


1 


2 H 


7 


(222) 


(40) 


9 


2 N 


7 


(222) 


(40) 


2 


2 I 


3 


(210) 


(20) 


1 


2 


5 


(221) 


(31) 





2 I 


5 


(221) 


(20) 


2 


2 Q 


7 


(222) 


(40) 
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No initial estimates were provided to the program. Non-physical situ- 
ations may then arise in the SCF process [6]. The output shows one such 
situation during the first iteration, but the program recovers and during 
later iterasions such warnings no longer appear. The ratio of the potential 
and kinetic energy is close to -2.0, as required by the virial theorem, showing 
good convergence for the wave function. 
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TEST RUN OUTPUT 



»hf 



HARTREE-FOCK . 96 



THE DIMENSIONS FOR THE CURRENT VERSION ARE: 
NWF= 20 N0=220 



START OF CASE 



Enter ATOM, TERM, Z 

Examples: 0,3P,8. or Oxygen, AV, 8. 
>Nd+,4F,60. 

List the CLOSED shells in the fields indicated (blank line if none) 

etc. 

> Is 2s 2p 3s 3p 3d 4s 4p 4d 5s 5p 

Enter electrons outside CLOSED shells (blank line if none) 

Example: 2s(l)2p(3) 
>4f (4)6s(l) 

There are 13 orbitals as follows: 

Is 2s 2p 3s 3p 3d 4s 4p 4d 5s 5p 4f 6s 

Orbitals to be varied: ALL/N0NE/=i (last i) /comma delimited list/H 
>all 

Default electron parameters ? (Y/N/H) 

>y 

Ambiguous l**n parent term: Enter term and Nr for f-subshells 
>3F3 

Default values for remaining parameters? (Y/N/H) 

>y 

WEAK ORTHOGONALIZATION DURING THE SCF CYCLE= T 

SCF CONVERGENCE TOLERANCE (FUNCTIONS) = 1.00D-08 

NUMBER OF POINTS IN THE MAXIMUM RANGE = 220 



ITERATION NUMBER 1 



SCF CONVERGENCE CRITERIA (SCFTOL*SQRT(Z*NWF) ) = 2.8D-07 



C( Is 


6s) 


= 


-0.03084 


V( 


Is 


6s) 


= 


-1565.84714 


EPS 


= 


.000020 


C( 2s 


6s) 


= 


-0 . 06406 


V( 


2s 


6s) 


= 


-274.28004 


EPS 


= 


.000234 


C( 3s 


6s) 


= 


-0.16604 


V( 


3s 


6s) 


= 


-72.19584 


EPS 


= 


.002300 


C( 4s 


6s) 


= 


. 02646 


V( 


4s 


6s) 


= 


-10.07054 


EPS 


=-0. 


.002627 


C( 5s 


6s) 


= 


0.12630 


V( 


5s 


6s) 


= 


-0.97428 


EPS 


=-0. 


.129631 


E( 6s 


Is) 


= 


0.00000 


E( 


Is 


6s) 


= 


0.00000 








E( 6s 


2s) 


= 


-0.00002 


E( 


2s 


6s) 


= 


-0.00001 








E( 6s 


3s) 


= 


-0.00025 


E( 


3s 


6s) 


= 


-0.00012 








E( 6s 


4s) 


= 


0.00017 


E( 


4s 


6s) 


= 


0.00008 








E( 6s 


5s) 


= 


0.01744 


E( 


5s 


6s) 


= 


0.00872 









EL 



ED 



AZ 



NORM 



DPM 
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Is 


3128.6805851 


921.2342946 


1.0536970 


4 


19D-02 








2s 


548.0200388 


304.9174534 


1.0620778 


6 


60D-02 








2p 


522.7712426 


5018.8746437 


1.3675876 


1 


81D-01 








3s 


144.2704835 


142.1548417 


5.7746458 


4 


17D-01 








3p 


140.9608657 


2469.2031309 


0.9319004 


4 


25D-01 








3d 


119.7093213 


8458.4713990 


0.7705958 


5 


91D-01 








4s 


43.7551400 


86.7032048 


0.3218645 


2 


57D+00 








4p 


39.3455261 


1332.3542693 


0.6952517 


2 


38D+00 








4d 


26.7211399 


4199.0102081 


0.4802196 


1 


82D+00 








4f 


7.5878908 


1305.5948122 


0.5721312 


9 


82D-01 








5s 


6.3617680 


27.3266729 


1.0694868 


8 


96D-01 








5p 


4.0723277 


429.0689970 


0.9602206 


7 


93D-01 




ED = 


0.029648 


ADJUSTED TO 


ALLOWED MINIMUM 


ENERGY 












6s 


0.4200626 


8.7147854 


0.2148108 


4 


65D-01 








4p 


14.0173622 


1022.4876672 


2.1997080 


3 


93D-01 








4d 


8.3655861 


3338.2783677 


1.6493227 


2 


70D-01 








4s 


20.1101311 


65.3845359 


8.0426513 


5 


70D-01 








4f 


1.9641284 


1955.2655220 


1.1526275 


1 


99D-01 








5p 


2.6937391 


310.8797108 


0.8665884 


3 


37D-01 


< 


Is 


2s>= 


2.4D-04 












< 


Is 


3s>= 


1 . 2D-03 












< 


2s 


3s>= 


1.0D-02 












< 


2p 


3p>= 


3.4D-03 












< 


Is 


4s >= 


5 . 3D-04 












< 


2s 


4s>= 


8.0D-03 












< 


3s 


4s>= 


6.0D-02 












< 


2p 


4p>= 


6.3D-03 












< 


3p 


4p>= 


7.4D-02 












< 


3d 


4d>= 


6.5D-02 












< 


Is 


5s>= 


1 . 4D-04 












< 


2s 


5s>= 


3.2D-03 












< 


3s 


5s>= 


2.4D-02 












< 


4s 


5s>= 


5.6D-02 












< 


2p 


5p>= 


4.2D-04 












< 


3p 


5p>= 


1.1D-03 












< 


4p 


5p>= 


-1.4D-01 












< 


Is 


6s>= 


-3.7D-05 












< 


2s 


6s>= 


8.0D-04 












< 


3s 


6s>= 


1.0D-02 












< 


4s 


6s>= 


5.6D-02 












< 


5s 


6s>= 


4.0D-01 













Iterations Omitted for Brevity . . . 
ITERATION NUMBER 8 



SCF CONVERGENCE CRITERIA (SCFTOL*SQRT(Z*NWF) ) = 3.6D-05 



C( Is 


6s) 


= 


0.00000 


V( Is 6s) 


= 


-1527.02399 


EPS = 0.000000 


C( 2s 


6s) 


= 


0.00000 


V( 2s 6s) 


= 


-243.19942 


EPS = 0.000000 


C( 3s 


6s) 


= 


0.00000 


V( 3s 6s) 


= 


-54.06456 


EPS = 0.000000 


C( 4s 


6s) 


= 


0.00000 


V( 4s 6s) 


= 


-11.16158 


EPS = 0.000000 


C( 5s 


6s) 


= 


0.00000 


V( 5s 6s) 


= 


-1.49103 


EPS = 0.000000 
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E( 6s Is) = 


O.OOOOl E( 


Is 6s) 





00001 








E( 6s 2s) = 


0.00004 E( 


2s 6s) 





00002 








E( 6s 3s) = 


-0.00348 E( 


3s 6s) 


-0 


00174 








E( 6s 4s) = 


-0.00227 E( 


4s 6s) 


-0 


00113 








E( 6s 5s) = 


0.02248 E( 


5s 6s) 





01124 










EL 


ED 




AZ 


NORM 




DPM 






Is 


3054.8943660 


921 


1231536 


1.0000000 


5 


06D- 


-09 




2s 


487.2430204 


301 


6764224 


1.0000000 


1 


19D- 


-07 




2p 


460.9251144 


4957 


8113013 


1.0000000 


1 


18D- 


-07 




3s 


109.0072008 


137 


6751878 


0.9999999 


4 


08D- 


-07 




3p 


97.3900098 


2374 


4492227 


0.9999999 


4 


38D- 


-07 




3d 


76.0249753 


7959 


4441188 


0.9999999 


3 


39D- 


-07 




4s 


23.2551364 


65 


1851169 


1.0000000 


6 


20D- 


-07 




4p 


18.6793842 


1101 


2565739 


1.0000000 


6 


83D- 


-07 




4d 


10.5291882 


3652 


0923901 


1.0000002 


4 


94D- 


-07 




4f 


1.4033321 


1922 


9003308 


1.0000011 


1 


71D- 


-06 




5s 


3.8125602 


25 


1175225 


1.0000001 


6 


86D- 


-07 




5p 


2.4248720 


387 


7133668 


1.0000000 


7 


49D- 


-07 




6s 


0.7066424 


7 


8763382 


0.9999988 


9 


91D- 


-07 


< Is 


2s>= 7.2D-10 
















< Is 


3s>= 4.8D-10 
















< 2s 


3s>=-1.6D-09 
















< 2p 


3p>=-6.9D-10 
















< Is 


4s>= 1.8D-10 
















< 2s 


4s>= 1.8D-09 
















< 3s 


4s>=-1.2D-08 
















< 2p 


4p>= 1.6D-09 
















< 3p 


4p>=-7 . 7D-09 
















< 3d 


4d>=-1.5D-08 
















< 3s 


5s>=-8.5D-09 
















< 4s 


5s>=-2.1D-07 
















< 2p 


5p>=-1.2D-10 
















< 3p 


5p>=-9 . 8D-09 
















< 4p 


5p>=-2 . 2D-07 
















< 2s 


6s>= 1.3D-10 
















< 3s 


6s>= 2.1D-10 
















< 4s 


6s>=-2.9D-08 
















< 5s 


6s>=-1.4D-07 
















TOTAL 


ENERGY (a.u.) 

















Non-Relativistic -9283.50837583 
Relativistic Shift -305.34173297 
Relativistic -9588.85010880 



Kinetic 9283.50834626 
Potential -18567.01672209 
Ratio -2.000000003 



Additional parameters ? (Y/N/H) 



Do you wish to continue along the sequence ? 



END OF CASE 
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